Coverage- adjusted confidence intervals for a 
binomial proportion 

Mans Thulin^ 



Abstract 

We consider the classic problem of interval estimation of a proportion p 
based on binomial sampling. The "exact" Clopper-Pearson confidence in- 
terval for p is known to be unnecessarily conservative. We propose coverage- 
adjustments of the Clopper-Pearson interval using prior and posterior dis- 
tributions of p. The adjusted intervals have improved coverage and are 
often shorter than competing intervals found in the literature. Using new 
heatmap-type plots for comparing confidence intervals, we find that the 
coverage-adjusted intervals are particularly suitable for p close to or 1. 
Keywords: Binomial distribution; Confidence interval; Proportion. 

1 Introduction 

Constructing a confidence interval for a proportion p based on a binomial sample 
is a basic but important problem in statistics. Due to the discreteness of the 
binomial distribution, it is not possible to construct confidence intervals with exact 
coverage. Thus an interval based on normal approximation, known as the Wald 
interval, is taught in virtually every introductory statistics course. The interval is 
p±Za/2\/f>4h^) where p = X/n is the sample proportion, q = 1—p and Za/2 is the 
100(1 — a/2)th percentile of the standard normal distribution. 

Numerous authors have remarked on the surprisingly poor performance of the 
Wald interval. Errors in the approximation due to discreteness and skewness (for 
small p) can have significant impact on the coverage of the interval even for large n. 
In recent years, its weaknesses have been thoroughly investigated in comparisons 
of confidence intervals for p. Brown et al. [31 H] gave examples of the erradic 
behaviour of the Wald interval, compared several intervals in terms of coverage 
and expected length and obtained general asymptotic results using Edgeworth 
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expansions. Pires & Amado [K] compared twenty methods using different criteria. 
For recent developments and discussions, see for instance [S], [3 [101 1121 ESI El]- 

A natural alternative to the Wald interval is the Clopper- Pearson interval |6]. 
It is based on the inversion of the equal-tailed binomial test and hence the interval 
contains all values of p that aren't rejected by the test at confidence level a. The 
lower limit is thus given by the value of pi such that 



E 



k=X 

and the upper limit is given by the pu such that 

X 



n—k 



a/2. (2) 



The computation of pl and pu is simplified by the following equality from Johnson 
et al. [15 : 

E (^) /(I - pT-' = r fit, X,n-X + l)dt. 



k=X 



where f{t,r,s) is the density function of a Beta{r,s) random variable. Conse- 
quently, the endpoints of the Clopper-Pearson interval Icp = {pl,Pu) are beta 
quantiles: 

I^p = (^B{a/2, X,n-X + 1), 5(1 - a/2, X + l,n-X) 

Icp is exact in the sense that the minimum coverage over all p is at least 1 — a. For 
most values of p however, especially values close to or 1, it is far too conservative, 
giving a coverage that is much larger than the nominal coverage. 

As several authors have pointed out [H [3], [13] it is often more natural to study 
the mean coverage rather than the minimum coverage. In this paper we construct 
coverage- adjusted Clopper-Pearson intervals with the mean coverage in mind, com- 
bining Bayesian and frequentist reasoning. The intervals are adjusted to have mean 
coverage 1 — a with respect to either a prior or a posterior distribution of p. The 
corrected intervals are seen to have several desirable properties in the frequentist 
setting. 

A class of coverage- adjusted Clopper-Pearson intervals is introduced in Section 
[2j In Section [3] these intervals are compared to other popular intervals and new 
heatmap-style plots for comparing confidence intervals are introduced. The text 
concludes with a discussion in Section |4] and an appendix with proofs, tables and 
several figures. 
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Figure 1: Actual coverage of the nominal 95 % Clopper-Pearson interval. 



2 Coverage-adjusted Clopper-Pearson intervals 
2.1 Definition 

As has already been mentioned, Iqp is often unnecessarily conservative. This is 
illustrated in Figure [TJ It is clear from the figure that if we are willing to accept an 
interval which has a coverage less than 1 — a for some values of p, the performance 
of Icp can be improved by choosing a larger a, in which case the actual coverage 
would be closer to the desired coverage for most values of p. The question, then, is 
how to choose the new a. We propose that a' should be chosen to satisfy a mean 
coverage criterion. 

Definition 1. Let f{-) be a density function on (0,1). A mean coverage corrected 
1 — a Clopper-Pearson interval Iqcp = {Pl,Pu) is given by the unique solution to 



I — n V / 



where a' satisifies 



C{a',n)= [ Pipelcp)- f{p)dp 
Jo 

= rEl(pe/cp(X,aO)f 

Jo v-n V 



(3) 

p^{l-py^-^.f(p)dp = l-a, 



x=o 

i.e. a' is such that the mean coverage of Igcp with respect to f is 1 — a. 
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Note that this simply is the ordinary 1 — a' Clopper-Pearson interval, with a' 
chosen so that the mean coverage is 1 — a. What differs is that a' needs to be 
determined before the endpoints are computed. 

It should be pointed out that the adjusted intervals inherit important prop- 
erties from Icp- They are fully boundary-respecting, so that Igcp ^ (0,1), and 
equivariant in the sense of [2] , meaning that the corresponding interval for 1 — p 
is (1 — pu, 1 — Pl)- Furthermore, they have very favourable location properties in 
terms of the Box-Cox index of symmetry and balance of mesial and distal non- 
coverage, as described by Newcombe [H]. Finally, the minimum coverage over all 
p is guaranteed to be at least 1 — a'. 

The choice of / affects the performance of Igcp greatly. / can be thought of as 
a weight function on (0, 1), used to put more weight on the performance for certain 
parts of the parameter space. In the following, we will refer to / as being either 
a prior or posterior density, to show the connection between this weight function 
and Bayesian ideas. 



2.2 Prior mean coverage corrections 

The use of a prior distribution / for coverage-adjustments can be motivated by 
the fact that in virtually all investigations, the experimenters will have some prior 
idea about how large p is. In particular, it is often clear beforehand if p is close to 
or far away from 1/2. 

Icp is symmetric in p in the sense that the interval has the same properties 
for p and 1 — p. For this reason, it is reasonable to use a symmetric prior for p. 
Beta{r, r) priors, being conjugate priors of the binomial distribution, are a natural 
choice here. We divide the parameter space into three cases: 

p close to or 1. When p is small, a prior with r < 1 should be used, 
as such priors put more weight on the tails of the distribution. We will use the 
Beta{l/2,l/2) prior in the following, but smaller r can certainly be used. The 
coverage-adjustments will generally be larger for small p, as the overcoverage of 
Icp is largest in this part of the parameter space. 

p close to 1/4 or 3/4- For medium-sized p, we wish to put approximately 
the same weight on the tails and the centre of the distribution. The uniform 
Beta{l, 1) prior is ideal for this. The resulting interval will however give a slight 
undercoverage for p closer to 1/2, so if there is some worry that that p may be 
above 0.40, say, a prior with r slightly greater than 1 could be used. The interval 
constructed using the uniform prior seems to coincide with a corrected interval 
that was described informally by Reiczigel pSl . 



p close to 1/2. If p is believed to be closer to 1/2, a prior with r > 1 is 
recommendable. We will use the Beta{2, 2) prior. The coverage-adjustments will 
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be smaller in this part of the parameter space, as Icp comes closest to attaining 
its nominal coverage around p = 1/2. 




Figure 2: Three useful prior distributions for p. 



2.3 Posterior mean coverage corrections 

Having used priors for coverage correction, it seems natural to consider using a 
posterior distribution of p for coverage-adjustments, in order to get closer to the 
nominal coverage in areas of the parameters space that given the data are more 
likely to contain p. 

With a Beta{r, s) prior for p, the posterior distribution is Beta{X + r, n — X + s) 
with density function 

^(p) = Wt W 0<p<l, 

p[X + r,n — X + s) 

where ■) is the beta function. Thus the posterior coverage corrected Clopper- 
Pearson interval Igcp is, given X, determined by the condition ([s]) with the func- 



5 



tion 



/•I " /„\ X+r-l(i _ \n-X+s-l 

In the comparison later in the text, we will use the Beta{l/2, 1/2), Beta{l, 1) and 
Beta{2, 2) priors, with the same reasoning as in the previous section. 

Conditioning the coverage-adjustments on the data may seem hazardous in a 
frequentist setting, but as we will demonstrate in Section [3| this approach leads 
to short confidence intervals with good coverage properties. 

2.4 Determining the adjusted confidence level 

While a' can be approximated by using an asymptotic expansion for the coverage 
to solve the equation (|3| approximately, it is more convenient to use a numerical 
method with exact coverages. The following lemma ensures that C is continuous 
and decreasing in a'. This guarantees that a' easily can be found numerically by 
using for instance bisection to solve the equation C{a', n,r, s) = 1 — a. The proof 
of the lemma is given in the Appendix, along with a table of a' for different choices 
of a and n for f{p) = 1. 

Lemma 1. Let Icp{X,a) = {pL{X,a),pu{X,a)) be the 1 — a Clopper-Pearson 
interval and let f{p,r,s), < p < 1, be the density of the Beta{r,s) distribution. 
The mean coverage of Icp{X,a) with respect to the density f{p), 

C{a, n, r, s) = j^' J^^Up e Icp{X, a)) (^^^''(l - pT^^'fiP, r, s)dp, 

is continuous and strictly decreasing in a. 

The algorithm for finding a' using bisection is as follows. 
Algorithm. Given a tolerance tol, a, n and a density /: 

1. Start with an initial lower bound a'^Q = a and an upper bound a'uQ. The 
initial guess is = {a'j^ q + a'jjQ)/2. 

2. Set i = 0. 

3. While |1 - a - C{a[,n)\ > tol: 

• If C{a'i, n)>l-a then a'^^i^^ = a[, a'jj^^^^ = a'^^^ and a^+i = [a'^^i+i + 
«W+i)/2- 
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• Else = a'^^^+i = ^ and a[_^^ = {a'^^i+i + a[7,,+i)/2. 

• i = i + 1. 

4. a' = a-. 

For the algorithm to converge, two conditions must be satisfied. First, a' < 
a^Q, i.e. a' must not exceed the upper bound. Second, C(a^, n) must be computed 
with sufficient precision {tol determines what is sufficient). 

Implementations of the above algorithm in R and MS Excel are available from 
the author. 



2.5 An example 

We illustrate the use of the coverage-adjustments with clinical data from an in- 
ffuenza vaccine study performed by Heinonen et al. [9]. n = 96 fully vaccinated 
children younger than 2 years were included in the study. X = 4 of these con- 
tracted inffuenza during the 2007-08 inffuenza season. 

The 95% Clopper-Pearson interval for the proportion of vaccinated children 
younger than 2 years that will contract inffuenza is (0.012,0.103). Using a prior 
Beta{l, 1) correction, we get a' ~ 0.06967. Letting B{-,r,s) be the quantile func- 
tion of the Beta{r, s) distribution, the coverage-adjusted conffdence interval is 

(^5(0.06967/2,4,93), 5(1 -0.06967/2,5,92)) = (0.013,0.098). 

Using a posterior Beta{l/2,l/2) correction, a' ~ 0.09385 and the interval is 
(0.014,0.094). 



3 Comparison of intervals 
3.1 Other intervals 

Following the comparison performed by Brown et al. [3J, two conffdence intervals 
for p have emerged as being the intervals to which all other intervals should be 
compared. These are the Wilson and Jeffreys prior intervals, presented next. 

The Wilson interval. Like the Wald interval, the Wilson [17J score interval is 
based on an inversion of the large sample normal test 



p — p 



d{p) 



< z, 



a/2, 



where d{p) is the standard error of p. Unlike the Wald interval, however, the 
inversion is obtained using the null standard error {pq/nY^'^ instead of the sample 
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standard error (pg/n)^/^. The solution of the resulting quadratic equation leads 
to the confidence interval 

J ^ + Zc,/2 [71 , 2 lA 

Iw = ± — -^^\ Pqn + z^,J4:. 

n + zli^ n + V "/^ 

ly/ typically has coverage close to the nominal coverage and comparatively short 
expected length. Indeed, it can be shown |1] that Iw has some near-optimal length 
properties among intervals with nominal coverage 1 — a. Iw is therefore the natural 
benchmark for new confidence intervals. 

The main drawback of Iw is that its coverages oscillates too much for fixed 
n and p close to or 1. Recently, Guan [8] proposed a small modification of 
the interval that solves this problem. The improved coverage comes at the cost 
of a slightly wider interval. As the results of our comparison would not change 
qualitatively if the Guan interval were to be used instead of Iw-, we stick to the 
more familiar unmodified version. 

The Jeffreys prior interval. Let X ~ Bin{n,p) and let p have prior distri- 
bution Beta{r, s). Then the posterior distribution is Beta{X + r,n — X + s) and 
letting B{a, r, s) denote the a-quantile of the Beta{r, s) distribution, a 100(1 — a)% 
Bayesian interval is 

(^B{a/2, X + r,n- X + s), B{1 - a/2, X + r, n - X + s)) . 

Pires & Amado [T3] used the uniform prior r = s = 1 in their comparison, whereas 
Brown et al. |3] used the Jeffreys prior r = s = 1/2. The difference between the 
two intervals is small. We use the latter and denote it Ij. 

I J has performance close to that of the Iw, and is often prefered when p is 
believed to be close to or 1. 

3.2 Numerical comparison 

In Figure [3]the actual coverages of some confidence intervals with nominal coverage 
95% are shown for p G (0, 0.5] when n = 25. All intervals are equivariant in the 
sense that the coverage is the same for p and 1 — p. Note that the coverages have 
been computed exactly (up to machine epsilon) and thus haven't been obtained 
by simulation. 

The Wilson interval has fairly good coverage properties when p isn't close to 0, 
in which case it oscillates wildly. The Jeffreys prior interval has similar coverage, 
except for one big dip for a moderately sized p. The coverage-adjusted intervals 
tend to have good coverage properties in the areas dictated by the prior distribution 
used for the adjustment. Some of the prior corrected intervals suffer from either 
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undercoverage or overcoverage for the parts of the parameter space that / put low 
weight on. 

The expected length of the intervals are shown in Figure |4j In many cases, 
the corrected intervals have shorter expected length than the Wilson and Jeffreys 
prior intervals. For some intervals, this is due to undercoverage caused by / putting 
more weight on a different part of the parameter space, but in some cases it is a 
consequence of a succesful correction. 

In order to compare intervals over the entire parameter space for different 
values of n, we use heatmap-type plots in Figures [5 10 Studying the plots for 



coverage and expected length at the same time gives us a good way of comparing 
the intervals. The heatmap-type plots, combined with more traditional plots such 
as those in Figures |3]j4| give a more complete comparison of different intervals than 
what has previously been possible. 

The Wilson interval is compared to the prior corrected Clopper-Pearson Beta{l, 1) 
interval in Figures |5]j6] and to the posterior corrected Clopper-Pearson Beta(0.5, 0.5) 
interval in Figures [^jS} The corrected intervals simultaneously offer greater cover- 
age and shorter expected length for small p. The difference is larger for small n 
and is particularly noticeable at the 99 % confidence level. 

In the comparison between the Jeffreys prior interval and the posterior cor- 
rected Clopper-Pearson Beta{0.5,0.5) interval in Figures |9 10, the corrected in- 
terval is found to offer at least as short intervals with the same actual coverage as 
the Jeffreys prior interval. 



4 Discussion 
4.1 Conclusions 

We introduced coverage- adjusted Clopper-Pearson intervals, where the intervals 
are adjusted to give mean coverage 1 — a with respect to either a prior or posterior 
distribution of p. We investigated the properties of several such intervals. The 
numerical results were presented graphically, partially with new heatmap-type 
plots. 

In the comparison with the benchmark Wilson and Jeffreys prior intervals, 
we found the coverage- adjusted Clopper-Pearson intervals to be preferable if p is 
believed to be close to or 1, as these intervals have both better coverage and 
shorter expected length in this setting. We have thus seen that it is possible to 
improve upon the Wilson and Jeffreys prior intervals for p close to or 1, if we 
are willing to accept that we use intervals that may have bad coverage properties 
in regions of the parameter space that are far from where our prior information 
indicates that p is. 
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In conclusion, the coverage- adjusted Clopper-Pearson intervals seem to be 
strong competitors against other methods for constructing confidence intervals 
for small binomial proportions. For p closer to 0.5, the Wilson interval seems to 
be preferable. 

4.2 Further developments 

The extension of the ideas presented here to one-sided intervals and to other distri- 
butions, such as the Poisson and negative binomial distributions is straightforward. 
Likewise, it should be possible to apply such corrections to tests about the differ- 
ence of two binomial proportions. It remains to be seen whether the corrections 
yield intervals with interesting properties in these cases as well. 

Apart from mean coverage, several other conditions can be used to ensure that 
the coverage is close to 1 — o; on average. Examples include median coverage con- 
ditions, minimum mean squared coverage error conditions and minimum absolute 
coverage error conditions. 
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A Appendix 



A.l Continuity and monotonicity of the mean coverage 



Since sums of continuous functions are continuous, it suffices to show that 



is continuous for fixed X, n, r and s. Seeing as f{p) = ^''^-'^(1 — pY^^/f3{r, s), this 
definite integral is a polynomial in Pl{,X, a) and pu{,X, a). Since the quantile func- 
tions of the beta distributions, and thus the limits of integration, are continuous 
in a, the continuity of C(a, n, r, s) in a follows. 

Similarly, as a) is strictly increasing in a and pu{X,a) is strictly de- 

creasing in a, and since — p)""'^ f{p,r, s) > for all p, the definite integral 
Q is strictly decreasing in a. C{a,n,r, s) is the sum of n + 1 strictly decreasing 
functions and thus also strictly decreasing. □ 



For a given prior distribution, a' is easily computed numerically given n, a and, in 
the case of a posterior correction, X. We give a table of a' for the prior corrected 
Clopper-Pearson Beta{l, 1) interval as an example. 



Proof of Lemma [I| Changing the order of summation of integration, the mean 
coverage can be rewritten as 





(4) 



A. 2 Tables 
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Table 1: ol for the prior corrected Clopper-Pearson Beta{l, 1) interval. 



n 


a = 05 


a = 0.01 


n 


a = 05 


a = 0.01 


Tl 


a = 05 


a = 0.01 


5 


0.1772 


0.0516 


55 


0.0769 


0.0171 


110 


0.0682 


0.0147 


10 


0.1280 


0.0331 


60 


0.0756 


0.0167 


120 


0.0674 


0.0145 


15 


0.1095 


0.0269 


65 


0.0745 


0.0164 


130 


0.0666 


0.0143 


20 


0.0995 


0.0237 


70 


0.0735 


0.0161 


140 


0.0660 


0.0141 


25 


0.0931 


0.0218 


75 


0.0726 


0.0159 


150 


0.0654 


0.0139 


30 


0.0885 


0.0204 


80 


0.0718 


0.0156 


160 


0.0649 


0.0138 


35 


0.0851 


0.0194 


85 


0.0710 


0.0154 


170 


0.0644 


0.0137 


40 


0.0825 


0.0186 


90 


0.0704 


0.0153 


180 


0.0640 


0.0135 


45 


0.0803 


0.0180 


95 


0.0698 


0.0151 


190 


0.0636 


0.0134 


50 


0.0785 


0.0175 


100 


0.0692 


0.0150 


200 


0.0632 


0.0133 



Next, we give some examples of the uncorrected, prior corrected Beta{l, 1) and 
posterior corrected Beta{l/2,l/2) Clopper-Pearson intervals for a few combina- 
tions of n, X and a. The posterior corrected Beta{l/2, 1/2) tends to get a larger 
correction, and thus shorter intervals, than the prior corrected Beta{l, 1) interval 
if the observed X is close to or n and a smaller correction if X is close to n/2. 
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Table 2: Some examples of coverage-adjusted Clopper-Pearson intervals. 





No correction 


Prior Beta{l, 1) 


Posterior 




1 b 




(T — 05 


— 01 


— 05 


— 01 


a = 0.05 


a = 0.01 


on 


1 
i 


nni o 
(^U.UUiz, 


(c\ nnnQ 
(^U.UUUo, 


(^U.UUZO, 


( r\ nnnfi 
(^U.UUUO, 


(0 0049 


(0 001 






0.2487) 


0.3171) 


0.2163) 


0.2815) 


u.iyzo J 


n Qi^Qi ^ 
u.zoyi ) 




O 
Z 


(c\ ni 




/"n ni fin 

\^U.UioU, 


(r\ nnfiQ 


/'n non7 


((\ nnoo 
(^u.uuyy, 






0.3170) 


0.3871) 


0.2829) 


0.3509) 


u.zDy/ ) 


U.OOD4 j 







(^U.UoOO, 


(^U.UOoo, 


/n 1 nQO 
(^U. iUoy, 


/n n7i s 
[v.v i io. 




(r\ n7nc; 
/ Uo, 






0.4910) 


0.5598) 


0.4559) 


0.5248) 


n /I fins') 
U.40Uo ) 


U.OZoi J 




1 n 
iU 


(^U.Z / zu, 


(^U.Zi I ( , 


( r\ Qni 7 
(^U.oUi ( , 


(r\ OA A 1 
(^U.Z44 ( , 


(^U.zyzy, 


/n oQfin 






0.7280) 


0.7823) 


0.6983) 


0.7553) 


n 7n7i 


n 7fi'? 1 


ou 


i 


(^U.UUUO, 




{r\ nnnfi 


l(\ nnno 

(^U.UUUZ, 


(0 001 6 


fO 0004 






0.1065) 


0.1394) 


0.0967) 


0.1282) 


n nfil o\ 

U.UoiZ ) 


nil 90^ 

U. i iZZ j 







(^U.Uooo, 


(^U.UZZ, 


( r\ nQ7R 
(^U.UO (0, 


(^U.UZOO, 


( ri nQS7 

(^U.Uc)o ( , 


(c\ nofi7 
(^U.UzD r , 






0.2181) 


0.2580) 


0.2058) 


0.2448) 


n ono7^ 
u.zuz / ) 


n o/ino^ 

U.Z4UZ J 




1 o 


(^U.ioUD, 


(^U.iUOD, 


/"n 1 QOc: 


/'nil QQ 


( n 1 Q7Q 
[U.iOlO, 


/nil on 
(^U.iizU, 






0.3817) 


0.4255) 


0.3676) 


0.4112) 


n Q7no^ 

U.o / UZ J 


n /1 1 Qfi^ 

U.4ioD ) 




oc; 
zo 




/'n Qi 


(^U.OOoO, 


(C\ QOfiO 
(^U.oZoZ, 


^U.oD44, 


(^U.oz4z, 






0.6447) 


0.6845) 


0.6314) 


0.6718) 


U.DoOD ) 


n fi7c;7^ 

U.D / / J 


100 


1 


(0.0003, 


(0.0001, 


(0.0004, 


(0.0001, 


(f) nnns 


fO 0002 






0.0545) 


0.0720) 


0.0508) 


0.0677) 


0.0413) 


0.0577) 
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(0.0490, 


(0.0382, 


(0.0518, 


(0.0405, 


(0.0523, 


(0.0410, 






0.1762) 


0.2020) 


0.1705) 


0.1959) 


0.1695) 


0.1946) 




25 


(0.1688, 


(0.1477, 


(0.1739, 


(0.1525, 


(0.1728, 


(0.1515, 






0.3466) 


0.3769) 


0.3396) 


0.3698) 


0.3410) 


0.3712) 




50 


(0.3983, 


(0.3689, 


(0.4052, 


(0.3756, 


(0.4031, 


(0.3735, 






0.6017) 


0.6311) 


0.4948) 


0.6244) 


0.5969) 


0.6265) 
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A. 3 Figures 



Prior correction: Beta(0.5,0.5) 



Prior correction: Beta(1,1) 



Prior coriection: Beta(2,2) 
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Posterior correction: Beta(0.5,0.5) 



Posterior correction: Beta(1,1) 
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Posterior correction: Beta(2,2) 
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Figure 3: Coverage of nominal 95 % intervals for n — 25. 




Figure 4: Expected length of nominal 95 % intervals for n — 25. 



Expected length of 95 % intervals Expected length of 99 % Intervals 




Figure 5: Comparison of the expected length of the Wilson and prior corrected 
Clopper-Pearson Beta{l, 1) intervals for different n and p. In the black points, 
the prior corrected interval has shorter expected length. In the grey points, the 
intervals have equal expected lengths. Based on a grid of 500 equidistant values 
of p in (0.001,0.5). 
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Figure 6: Comparison of the coverage of the Wilson and prior corrected Clopper- 
Pearson Beta{l, 1) intervals for different n and p. In the black points the prior 
corrected interval has greater coverage, in the white points the Wilson interval has 
greater coverage and in the grey points the intervals have equal coverage (when 
rounded to 3 decimal places). 
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Figure 7: Comparison of the expected length of the Wilson and posterior corrected 
Clopper-Pearson Beta{(}.5. 0.5) interval intervals for different n and p. In the black 
points, the posterior corrected interval has shorter expected length. In the grey 
points, the intervals have equal expected lengths. 
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Figure 8: Comparison of the coverage of the Wilson and posterior corrected 
Clopper-Pearson Beta{0.5, 0.5) intervals for different n and p. In the black points 
the posterior corrected interval has greater coverage, in the white points the Wil- 
son interval has greater coverage and in the grey points the intervals have equal 
coverage (when rounded to 3 decimal places). 
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Figure 9: Comparison of the expected length of the Bayesian Jeffreys prior and 
posterior corrected Clopper-Pearson Beta{0.5,0.5) intervals for different n and p. 
In the black points, the posterior corrected interval has shorter expected length. 
In the grey points, the intervals have equal expected lengths. 
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Figure 10: Comparison of the coverage of the Bayesian Jeffreys prior and posterior 
corrected Clopper-Pearson Beta{0.5, 0.5) intervals for different n and p. In the 
black points the posterior corrected interval has greater coverage, in the white 
points the Bayesian Jeffreys prior interval has greater coverage and in the grey 
points the intervals have equal coverage (when rounded to 3 decimal places). 
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